Astronomy & Astrophysics manuscript no. mennekens 


© ESO 2008 


May 6, 2008 





On the origin and evolutionary state of RZ Cas, KO Aql and S Equ 

(Research Note) 

N. Mennekens^, J. -P. De Greve^ W. Van Rensbergen^, and L. R. Yungelson^ 

' Astrophysical Institute, Vrije Universiteit Brussel, Pleinlaan 2, 1050 Brussels, Belgium 

e-mail: nmenneke@vub.ac.be, jpdgreve@vub.ac.be, wvanrens@vub.ac.be 
^ Institute of Astronomy of the Russian Academy of Sciences, 48 Pyatnitskaya St., 1 19017 Moscow, Russia 

e-mail: lry@inasan.ru 

Received February 7, 2008; accepted April 30, 2008 

ABSTRACT 

Aims. Determination of the present evolutionary state and the restrictions on the initial mass ratios of RZ Cas, KO Aql and S Equ. 
Methods. Comparison of mass gaining stars with evolutionary models of single stars with the same mass and subsequent comparison 
with accretion tracks from simultaneous conservative binary evolution. 

Results. The gainers are in an early main sequence stage (X^ > 0.5), with KO Aql being almost unevolved (assuming quasi-thermal 
equilibrium). The initial donor/gainer mass ratios Mji/Mgi must be larger than three to obtain the present mass and luminosity of the 
gainers. 

Key words, stars: binaries: eclipsing - stars: binaries: close - stars: individual: RZ Cas, KO Aql, S Equ 



1. Introduction 

RZ Cas, KO Aql and S Equ are semi-detached eclipsing Algol- 
type binaries for which we investigate restrictions on their ori- 
gin coming from the comparison of observations with theoretical 
evolution models. For RZ Cas we use the orbital and physical pa- 
rameters compiled by Maxted et al. (1994)| and the M given by 
[Lehmann & Mkrtichian (2007)| The M for RZ Cas was derived 
from O-C curves. For KO Aql and S Equ, we use the results 
of Soydugan et al. (2007)[ All these are rendered in Table 1 . To 
avoid any confusion, we label the hotter, mass-gaining compo- 
nent (commonly called the primary) as gainer, while the cooler, 
mass-losing component (the secondary) is labeled as donor. 
The mentioned authors conclude from their analysis that the 
three systems are showing prominent features of mass transfer. 
Together they form an interesting set from an evolutionary point 
of view. The total masses range from 2.94Mo, over 3.08Mo to 
3.66Mq, the mass ratios from 0.13 over 0.22 to 0.33. The total 
mass, gainer mass and period of RZ Cas are the smallest, but 
combine with the largest mass ratio, whereas this is the oppo- 
site for S Equ. KO Aql lies in the middle of the two. This might 
indicate that the three systems are in subsequent stages of mass 
transfer. The low mass ratios point towards an advanced state of 
mass transfer. 



2. Comparison with single star evolution 

A comparison of the location of the gainers with single star 
evolutionary tracks can teach us how far they deviate from 
normal stars. Single star evolution (initial composition: X - 
0.7, Z = 0.02) was computed for the three masses of the 
gainers in Table 1 with the same evolutionary code used by 



De Loore & Van Rensbergen (2005) The results are shown in 



Figure 1, together with the positions of the three mass gaining 



stars. The mass of the lowest track corresponds to that of RZ Cas, 
the next to that of KO Aql, and the highest to the mass of S Equ. 

It is clear that the observed gainers are underluminous com- 
pared to the theoretical counterparts, as was previously noted by 
[Soydugan et al. (2007)| The evolution of a close binary can be 
divided into two distinct phases: an initial phase of fast mass 
transfer (increasing M) on a thermal timescale, and a later phase 
of slower mass transfer (decreasing M) on a nuclear timescale 
(e.g. advanced state with well reversed mass ratio). 

Given the well reversed mass ratios (RZ Cas: 0.33; KO Aql: 
0.22; S Equ: 0.13) and the relatively low mass loss rates deduced 
from the period changes (~ 10"^ to 10^'^ M^/yr), we assume that 
the systems are well beyond the peak phase of mass transfer. In 
that case, the gainer is close to thermal equilibrium and the evo- 
lution is mostly determined by the inner structure. Thus, its po- 
sition in the main sequence band reflects the amount of nuclear 
burning that has taken place and we can derive a rough estimate 
for the amount of hydrogen by mass in the core (Xc). We find 
that for RZ Cas X^ ~ 0.6, for KO Aql X, ^ 0.7, and for S Equ 
Xc ~ 0.5. In conclusion, this points towards a ZAMS structure 
for KO Aql, and a more evolved structure for S Equ. 

3. Conservative considerations 

To explore the values for the masses and the period of the pro- 
genitor systems, we first investigate a possible conservative his- 
tory of mass transfer (conservation of mass and angular momen- 
tum). We specifically discuss the accretion behavior in some de- 
tail for the middle system KO Aql (in mass, mass ratio and pe- 
riod), and we will extend the conclusions to the other two. In the 
code, the structure of the gainer is calculated assuming that the 
added matter has the same entropy as the gainer and neglecting 
all effects of interaction. 

We calculated the system parameters of KO Aql backwards 
to initial states corresponding to qi = 1.01 (initial gainer mass 
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Table 1. Orbital and physical parameters of RZ Cas, KO Aql and S Equ 



Parameter 


RZCas 


KO Aql 


S Equ 




gainer donor 


gainer donor 


gainer donor 


P(d) 


1.195 


2.864 


3.436 


M{M^) 


2.21 ± 0.08 0.73 ± 0.02 


2.53 ± 0.05 0.55 ± 0.01 


3.24 ± 0.03 0.42 ±0.01 


R(Re) 


1.67 ± 0.03 1.94 ±0.03 


1.74 ±0.07 3.34 ±0.07 


2.74 ±0.09 3.24 ±0.10 


log T (K) 


3.934 3.672 


3.996 3.646 


4.049 3.721 


log L/Lo 


1.12 ±0.02 0.16 ±0.08 


1.41 ±0.06 0.56 ±0.06 


2.02 ±0.10 0.86 ±0.10 


M{MJyr) 


~ 1.3 • 10-« 


~ 2.6 • 10-'' 


~ 4.0 • 10"* 



2.5 



2.3 - 



2.1 



® 1 9 



o 1.7 - 



1.5 



1.3 



1.1 



4.12 




3.84 



Fig. 1. HRD with the positions of the gainers of S Equ, KO Aql, RZ Cas, and tracks of single star evolution for the mass values of 
Table 1. Central hydrogen abundance by mass, as interpolated from the models, is indicated next to the gainers (see text). 



Mgi = 1.53Mo), 1.32 (Mgi = 1.32Mo), and 1.61 (Mgi = 
1.1 8Mq), subsequently followed by a computation of the evo- 
lution of the gainers (through simultaneous evolution of both 
components). All three failed to reproduce the current param- 
eters of the gainer of KO Aql, in mass as well as in HRD posi- 
tion. The calculations were stopped when the mass transfer rate 
fell below 10 "Mo/yr, and nuclear evolution took over (and the 
donor returned to thermal equilibrium). The mass of the gainer 
at that moment for all three systems was in the range I.SOMq 
to I.86M0. The results for the gainer in the HRD are shown in 
Figure 2. 

All three resulted in gainers with masses and HRD positions 
far below the observed one. We note that the larger the initial 
mass ratio, the larger is the mass of the gainer and its luminos- 
ity in the slow phase. Therefore, we calculated a new set of two 
with more extreme mass ratios (3.40, Mgj = O.TOMq and 3.46, 
Mgi = O.69M0). The results are also shown in Figure 2. Both 
tracks show a similar behavior. As a result of the large differ- 
ence in thermal equilibrium timescales between the components 
(tKH,g/tKH,d ~ 3-5), the star first strongly deviates from thermal 
equilibrium, and strongly increases in volume. After mass ratio 



reversal the mass transfer rate decreases and the gainer evolves at 
nearly constant luminosity towards the ZAMS. Near the ZAMS 
the star resumes quasi-thermal equilibrium and moves down- 
ward along the ZAMS, finally setthng itself in the slow phase 
of mass transfer on a track corresponding to its mass and in- 
ner structure. The mass transfer rate evolves from a fast phase 
of ~ 1.7 ■ IQ^^yr with a maximum of a few times lO^^'Mo/yr 
to an average of IQ^^^M^/yr in the rest of the main sequence 
phase. The first series {Mgi = O.TOMq) results in an accretion 
track that runs through the observed position of KO Aql in the 
HRD. At that position it does not result in a mass Mg match- 
ing the observed one. Indeed, at the moment that it reaches the 
HRD position of KO Aql, its mass is 2.25Mo. During the further 
evolution, the mass transfer (of about M ~ 10"^ - 10"** Mg/yr) 
increases that mass to ~ 2.54M0. At M = 1.53Mg, it is already 
well evolved into the main sequence, its location corresponding 
to single star models with the same X^. 

The second track (Mg, = 0.69Mo) hes well above the ob- 
served position, but the higher mass loss rate results in the mass 
of KO Aql and even beyond. The characteristics of the model 
at the moment that Mg = 2.53Mq are log T^ff = 4.035 and 
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Fig. 2. HRD with accretion tracks for the gainer of KO Aql. The initial mass is indicated at the beginning of the track. Crosses 
indicate the positions where the masses and period of the systems equal those of KO Aql. 



log L/Lq - 1 .618. For the two calculations that attain the masses 
and period of KO Aql, the calculated properties at that time are 
given in Table 2. 

The different behavior of the two groups of systems can be 
understood in terms of the accretion timescale Mg/M versus the 
thermal timescale Ikh (~ M^/RL). This ratio is always larger 
than one for the mass ratios between one and two (at maximum 
mass loss rate it reaches ~ 6). However, for large initial mass 
ratios, the ratio of the timescales becomes as small as 0.01 at 
maximum M, while it is ~ 100 when Mg = M„hserved- 

In the slow phase of mass transfer, all calculated series have 
gainers that are close to thermal equilibrium. For the large ini- 
tial mass ratios, the ratio L„ucl Lsurjace ~ 0.98 when the mass 
of the gainer matches the observed one (compared to a value 
of 0.05 during maximum mass loss). For the three systems with 
mass ratios between one and two we obtain average values of 
0.91 (slow phase) and 0. 16 (maximum mass loss rate). When ac- 
cretion starts, since the thermal timescale of the donor is much 
smaller than for the gainer, the core is compressed almost adia- 
batically. As a result the temperature and density rise and gen- 
erate an energy surplus in the interiors. This results in an in- 
crease of L and R. The difference between the stars is in the abil- 
ity of their external layers to absorb this flux of energy. Similar 
results were found by ' Webbink ( 1 976)| and [Webbink ( 1 977)1 I" 
more massive stars, thermal timescales are similar and the gainer 
evolves more or less along the TAMS, while in lower mass stars 
it deflects from the ZAMS. 

Hence, to obtain an accretion model near the observed HRD 
position of the gainer, an initial mass ratio of ~ 3.4 is required. 
However, as the figure shows, during the fast mass transfer the 
gainer moves almost vertically in the HRD, obtaining a radius 
that largely exceeds the Roche radius. The conservative evolu- 



tion beyond the point where Rg - Rcit is therefore formal. Thus, 
a non-conservative (liberal) approach is necessary. 

The other two systems give similar conclusions. To reach 
the observed luminosity of the observed gainer, one has to start 
from a large mass ratio. The position of RZ Cas is reached 
with qi - 2.07 (Mg, = O.7O5M0), but then the model mass 
is only 1.89Mo. The observed mass is reached for qi - 2.94 
{Mgi = 0.705Ma), but with log T,ff = 3.98 and log L/Lq = 1.27. 
Here also, a large radius excess over the Roche radius occurs 
during the first and fast phase of mass transfer. Similarly, for 
S Equ a value of ^, ~ 3.3 - 3.4 is minimally needed, with even 
greater radius excess. For S Equ an additional problem occurs 
with conservative evolution. The large mass ratio results in large 
values for M^i, with radii that are already matching the Roche 
radii at the beginning of the main sequence evolution. 

The observed underluminosity remains a problem. The three 
systems with mass ratio between one and two were stopped at 
mass Mg - I.8OM0 to I.86M0. They turn out to be sHghtly 
underluminous. Taking for example the middle system, Mg, = 
1.32Mo, we obtain at Mg = 1.84Mo that log Lg/Lo = i.09 
{Xc - 0.54) whereas a 1.84Mq single star has log L/Lg = 1.12 
at the same Xc- However, in the large mass ratio series we find 
log L/Lq ~ 1.6 - 1.7 at Mg - 2.53M0 (equal to that of a single 
star of 2.53Mo), whereas KO Aql has log L/Lq = 1.4. 

4. Conclusion 

We have found that in order to reproduce the gainers of RZ Cas, 
KO Aql and S Equ through conservative binary evolution, one 
would have to start from large initial mass ratios (typically > 3). 
The corresponding high mass loss rates in the beginning of the 
RLOF result in radii larger than the Roche radii. Hence, a non- 
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Table 2. Parameters of the different calculations when they have the masses and period of KO Aql. 





0.69Mo 


O.VOMo 


log T, (K) 


4.035 


3.997 


log Lj/Lq 


1.618 


1.727 




8.1 • 10-* 


1.8- 10"' 


^nuCygl ^surface,g 




0.98 




0.70 


0.48 



conservative approach is needed, in which processes are con- 
sidered that enhance the period, without losing too much mass 
(as this would result in a higher initial mass ratio and problems 
in reaching the final mass and luminosity in the accretion pro- 
cess). We have planned to undertake such a study in the coming 
months, exploring also the reasons for the observed underlumi- 
nosity of the gainers. 
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